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Abstract 

A general algebraic condition for the functional independence of 2n — 1 con- 
stants of motion of an n-dimensional maximal superintegrable Hamiltonian 
system has been proved for an arbitrary finite n. This makes it possible to con- 
struct, in a well-defined generic way, a normalized Nambu bracket which pro- 
duces the correct Hamiltonian time evolution. Existence and explicit forms of 
pairwise compatible multi-Hamiltonian structures for any maximal superinte- 
grable system have been established. The Calogero-Moser system, motion of a 
charged particle in a uniform perpendicular magnetic field and Smorodinsky- 
Winternitz potentials are considered as illustrative applications and their sym- 
metry algebras as well as their Nambu formulations and alternative Poisson 
structures are presented. 
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I. INTRODUCTION 



Nambu mechanics is a generalization of the Hamihonian formulation of classical mechan- 
ics in that it replaces the usual binary Poisson bracket (PB) to higher order ?7,-ary bracket, 
generically called Nambu bracket (NB), and specifies the dynamics in terms of n — 1 "gen- 
eralized Hamiltonian" functions . The original motivation of Nambu was to show that 
the Hamiltonian mechanics is not the only formulation that makes a statistical mechanics 
possible. Relevance of Nambu mechanics to membrane theory has been put forward and 
a form of quantized Nambu mechanics has been purposed as a nonlinear generalization of 
geometric formulation of quantum mechanics [0,0. Nambu formulation may also give some 
insights into the theory of higher order algebraic structures and their possible physical sig- 
nificance [^]. Unfortunately, up to now only few examples of dynamical systems which admit 
Nambu formulation have been given. The Euler equations for three dimensional (3D) rigid 
body were the only example given by Nambu. Then, the equations of Nahm system in 
the theory of static SU{2) monopoles were realized in this formulation 0,^. Connection 
between Nambu mechanics and so-called superintegrable systems has been the subject of 
some recent studies [[7|,p|,p|,p!0| . 

A Hamiltonian system of n degrees of freedom is called to be completely integrable, in 
the Liouville-Arnold sense, if it admits n functionally independent, globally defined con- 
stants of motion in involution (i.e., commuting with respect to PB) |[T2| , p!3[| . A completely 
integrable system is called superintegrable if it allows k additional constants of motion. Not 
all constants of motion of a superintegrable system can be in involution but they must be 
functionally independent, otherwise the extra invariants are trivial. Superintegrability is 
said to be minimal if k = 1, and maximal if k = n — 1 |]l^|15lll6i . An nD maximally 



superintegrable Hamiltonian system can be specified by the following set 

SlH{2n - 1) = {H, H„ A, : {H, H,} = = {H„ H^}, {H, A,} = 0, F ^ 0} , (1) 

where H is the Hamiltonian of the system. Hi, Aj; i, j = 1, 2, . . . , n — 1 are the additional 
constants of motion and F denotes the following {2n — l)-form 

T = dH A dHi A ■ ■ ■ A dHn-i AdAiA--- A dA^-i- (2) 

Here d and A denote the usual exterior derivative and exterior product of Cartan calculus. 
For functional independence F must be different from zero on a dense subset of the underlying 
symplectic manifold endowed with the PB {, }. 

In this paper three main points concerning the fundamental structure of a Slni^n — 1) 
system for any finite n are established. We shall first prove that the constants of motion 
of a Slni^n — 1) system are functionally independent where (n — 1) x [n — 1) matrix B 
with elements Bij = {Hi, Aj} is nonsingular. Secondly we shall construct the normalized NB 
which produces the correct Hamiltonian time evolutions. This means that all nD maximally 
superintegrable Hamiltonian systems admit Nambu formulation. We then show that every 
Slni^n — 1) system admits 2n — 1 multi-Hamiltonian structures. Statements that we have 
proved make these facts possible, in a well defined generic way and independently from the 
forms of Hamiltonians. Multi-Hamiltonian structures of maximally superintegrable systems 
were considered for the first time in Ref. P] by using a different approach from ours. In this 
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context, our geometric proofs of Jacobi identity and of compatibility condition for alternative 
Poisson structures are different but validate and, in a sense, are complementary to that of 
Ref. i. 

In the next section main points of the Nambu mechanics, mainly those needed for the 
subsequent investigation are briefly reviewed ( for more details we refer to 0). In Sec- 
tion III a coordinate-free form of canonical NB and normalized NB are introduced, and 
the Jacobi identity for a normalized binary bracket induced from NB is established. The 
general algebraic condition for the functional independence is proved in section IV. Nambu 
formulation of a Slni'^n — 1) system is established in section V where we also point out 
some general facts concerning the structure of symmetry algebras of maximally superinte- 
grable systems. Multi-Hamiltonian structures are taken up in section VI. As applications the 
Calogero-Moser system, motion of a charged particle in a uniform perpendicular magnetic 
field (this will be referred to as (classical) Landau problem) and Smorodinsky-Winternitz 
potentials are considered in the final section where their symmetry algebras and explicit 
forms of their alternative Hamiltonian structures are established. 

Our further notational conventions are as follows. We shall denote the linear spaces of all 
vector fields and (differential) p-forms on a smooth manifold M of dimension n, respectively, 
by X{M) and Ap(M), < p < n. Together with their commutative (and associative) algebra 
structure with respect to usual point-wise product, the linear space of all smooth functions 
(0-forms) defined on M will be represented by A. We shall denote the vector fields by bold 
face letters, adopt the Einstein summation convention over repeated pair of contravariant 
and covariant indices and use the shorthand A^" = A ® ■ ■ ■ ^ A {n times). 

II. NAMBU MECHANICS 

NB of order n is the real multilinear map {, } : — > A which has, for all 
fj^Qj G A, the following properties. 

i. Skew-symmetry 

{/l, /n} = (-1)^{/<7{1), /a(n)}, (3) 

where a is a member of the permutation group Sn and e is the parity of permutation a 
{e = 1 for odd permutations, and e = for even permutations). 

ii. Derivation (the Leibniz rule) 

{/1/25/3! /n+l} = /l{/25 /s! /n+l} + /2{/l5 /s! /n+l}- (4) 

iii. Fundamental identity (a kind of generalized Jacobi identity) 

n 

{/l, • • • , /n-1, [gi, • • • , 9n}} = ^{gi, • • • , {/l, • • • , fn-1, Qk}, • • • , Qn} ■ (5) 

k=l 

With respect to NB, A acquires another algebra structure henceforth denoted by A^. 

Nambu dynamics is determined by n — 1 Hamiltonian functions hi, . . . , G An and 
is described, for any / G An, by the Nambu-Hamilton (NH) equations of motion 
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-77 — ^Nnif) — {/, hi, . . ., hn-l}, 



(6) 



where Xj^h is called the NH vector field corresponding to hi, . . . , /in-i- 

NB of order n induces infinite family of lower order NB, including the family of Poisson 
structures, all of which satisfy corresponding fundamental identities (FIs) that follow from 
(5). Below we shall concentrate only on the induced Poisson structures. For a fixed set of 
n — 2 Hamiltonian functions G we define the Nambu induced PB as follows 

{/, 9}np = {/, gji, • • • , /n-2}, (7) 
where /, g E An are arbitrary. If in Eq. (5) we take 

fn-l = f, ■ 9n-l = g, Qn = h, Qi = ff, i = 1, . . . , n - 2 

then, by virtue of (3), the first n — 2 terms at the right hand side (5) vanish and we get 

{/, {9, h}Np}NP + cp = 0, (8) 

where cp stands for cyclic permutations. Eq. (8) reveals the fact that {,}np satisfies the 
Jacobi identity. Note that all the fixed functions in the definition (7) are Casimirs of {, }np, 
that is, {fi, /Iatp = for all i = I, ■ ■ ■ , n — 2 and / G An- 



III. CANONICAL NB AND NORMALIZED NB 

The problem of constructing concrete realizations of NB is of great importance. In the 
case of M = R" the following form, called the canonical NB 

was provided by Y. Nambu. Here {x^, . . . , x") denote the local coordinates of R" and the 
right hand side stands for the Jacobian of the mapping / = (/i, . . . , /„) : R" — R". 

We shall now introduce a coordinate-free expression of the canonical NB that provides 
a considerable ease in proving the technical points of the paper. For this purpose we first 
associate the n — 1-form 

7 = (i/ii A ■ ■ ■ A dhn-i , (10) 

to 72 — 1 Hamiltonian functions hj G An- We then recall the Hodge map * : A^(R") — >■ 
A"~P(R") defined for any p-form w = {l/p\)wi^ ... i^dx^^ A ■ ■ ■ A dx^^ as follows fl^ 

*w = ——^ — -e*' ■■■ i Wi, i dx'^+^ A ■ ■ • A dx*", 

where Wi^,,,i^ are antisymmetric cotravariant components of w and e'^'' *" with e^---^ = 
1 is the nD completely antisymmetric Levi-Civita symbol. Note that with respect to A 
the Hodge map is linear and exterior product is bilinear. It is now obvious that, in local 
coordinates 
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^{df A 7) = {/, /.I, . . . , K-^} = m 

In that case, multilinearity, antisymmetry and derivation properties of the canonical NB 
are direct results of the linearity of the Hodge map with respect to A, and of the 
well-known properties of A-product and d. We can also associate the n — 1-form (3 = 
Ai ... in-i^x''^ A ■ ■ ■ A dx^"-^/{n — 1)! to the vector field f3 = {(3^, . . . , with components 
^ gfcji ...t„-ip,^ — 1)!. Then, (11) can also be written as 

*(rf/A7) = 7-V/, (12) 

where V stands for n-dimensional gradient operator and "■" denotes the usual inner product 
of R". The fundamental identity can be verified by taking / = {^fi, . . . , gn} in Eq. (12). 
We should note that Eq. (12) implies 

X^H = (-l)"-'^(7Arf) = 7-V, (13) 

for the NH vector field corresponding to hi, hn-i- As illustrative examples, let us 

consider the cases n = 2, 3. The vector fields corresponding to 'y = dh in the case of n = 2 
and 7' = dhi A dh2 in the case of n = 3 are easily found to be 

7 = id2h, -dih), 7' = VhxVh2 , 

where dj = d/dxj and X denotes the cross product of R^. The associated NH vector fields 
and canonical NBs can be written as follows 

XnhW) = d2hdj - dihd2f, X'^hU) = (V/iiX V/i2) ■ V/. 

The first is the usual PB of R^ and the second is the original NB first appeared in [|I|. In 
the next two sections we shall generalize these expressions for a Slni'^n — 1) system in the 
case of arbitrary n. 

To be precise, from now on we shall adopt, in accordance with related literature, the 
following definition : If the time evolution equations of a dynamical system can be written 
in terms of (canonical) NB then the system will be called to admit equivalent Nambu formu- 
lation. As it will be apparent in the next two sections in order to get the correct dynamics 
the induced PBs must be properly normalized. For this purpose, in terms of n — 2-form 

V = dfiA ■■■ A (14) 

we define 

{f,gyNP = C{f,gJi, fn-2} = C*idfAdgAr^), (15) 

where C G ^ is, for the time being, an arbitrary function. C will be referred to as the 
normalization coefficient and will be specified from the requirement that the NH equation 
produces the correct time evolution for any function. The generic form of C will be deter- 
mined in the next section, but before that it must be emphasized that {, }'^p satisfies the 
Jacobi identity for any C. To prove this let us consider 
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{h, {/, gYNpYNP = C'idh A d[C*idf AdgA r/)] A r]} 



(16) 



and its cyclic permutations for three functions h, f and g. Since in a nD space we have at 
most n functional independent functions and since in (15) we have, apart from C, n + 1 
functions, in the most general case one of them, say /, must be functional dependent to 
others. Hence we can take 

df = adh + bdg + ^ Cj/j, (17) 
1 

where a,b and Cj are arbitrary constants. On substituting this in (16) and in its cp we 
immediately see that their sum vanishes. 

Obviously, in the case of n = 2 we have C = 1. For known 3D examples, namely, for 
free rigid body and the Nahm equations we also have C = 1. However the requirement of 
nontrivial C is inevitable at least when n is an even integer greater than two. Although this 
normalization requirement have appeared in the literature, its generic form and important 
implications were not recognized. 



IV. AN ALGEBRAIC EXPRESSION FOR FUNCTIONAL INDEPENDENCE 

The phase-space of a Hamiltonian system is a 2n dimensional symplectic manifold M on 
which a symplectic structure is defined by a closed {dQ = 0) and nondegenerate symplectic 2- 
form Q. Two immediate implications of nondegeneracy are that [|11|,0]; (i) M is orientable 
with nowhere vanishing Liouville measure (volume form) 

Vl = 1^", (18) 

where Q'^ = Q A ■ ■ ■ A Q [n times) . (ii) There is a natural isomorphism between the vector 
fields and 1-forms defined by ^ ^ /i^ = i^Q, where : A^(M) —>■ Ap~^(M) is called the 
interior product operator defined for any p-form a by 

{i^a){^i,...,ip_i) = a(^,|i,...,|p-i). (19) 

If /ig is exact, that is, if fi^ = df, then ^ is called a Hamiltonian vector field corresponding to 
f & A and henceforth denoted by ^/ : i^j^Q = df. PB on M is defined by ^{$,f, $,g) = {/, g}- 
According to Darboux theorem at each point of M there are local canonical coordinates 
(g^, . . . , g", Pi, . . . , pn) in which Q takes the form Q = dq^ A dpj and leads us to the following 
coordinate expressions 

{/, g} = ^(O. ^g) = dgifdp^g - dpjdqjg, (20a) 
^f = dpj-dq,-dq,fdp^, (20b) 
VL = dq^ A--- Adq" AdpiA--- A dp^ . (20c) 

By definition of exterior forms and their A-products, the value of a A /5; a G A^, j3 E A'' 
on p + q vectors £ '^{M); k = 1, 2, . . . , p + q < 2n, is given, in the notation of |l^, by 
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(« A • • • , U,) = E • • • , ^^p)mn, ^.J> (21) 

Sp+q 

where ii < ■ ■ ■ < ip and ji < ■ ■ ■ < jg such that {ii, . . . , ip, ji, . . . , jg) is a permutation 
of (1, 2, . . . , p + g). The summation in Eq. (21) is over all permutations of the permutation 
group Sp+g, provided that the indices are partitioned into two ordered sets as given above. 

By making use of Eq. (21) we now evaluate the Liouville form Vl given by (20c) on 2n 
Hamiltonian vector fields A = 1, . . . , 2n, a.s follows 

= (— 1)"£*^ ^'^"dp^fi^ ... dp^Ji^dgifi^^^ ... dqnfi^^ 

/ -|^\n^(/"+l' • • • 5 /2n; /l; • • • ; fn) 

d{q\ g", Pi,..., p„) 

By the antisymmetry properties of determinant, or in view of Vl(^i, . . . , ^2n) = 
(— l)"VL(^n+i, • • • , ^2n, ^1, ■ • ■ , $n) , the abovc relation can be written as 

We now evaluate on the same set of Hamiltonian vector fields as follows 

= E(-imi,4} ■■■ {fk^jk,}, (23) 

where h < 12] ■ ■ ■ ;ki < k2 such that {ii, ^2, • • • , ^1, ^2) is a permutation of (1, 2, . . . , 2n), 
and we have used Eq. (20a). Let us make the identifications 

/i = /, /2 = if , f2+i = , fn+i+i = A; ^ = 1, 2, . . . , n - 1 , (24) 

where H, Hi, Ai G Slui^n — l) and / G .4 is an arbitrary function. In that case, at the right 
hand side of Eq. (23) only the terms in which H is paired with / and each Hi is paired with 
one of Aj give non-zero contributions. All other possible pairing are zero by very definition 
of maximal superintegrability. With this in mind let us consider a fixed partition 

{/il ) /i2 }{/?! 5 /j2 } • • • {/fel 5 /fc2 }• 

There are n different places for {/j^, fi^}, and for a fixed place of {/jj, fi^} there are n — 1 
possible places for {fj^, fj^}, and n — 2 places for the next pair (provided that the places of 



^When we were about to submit this paper we came across to very recent study [|T^ which derive 
Eq. (23) (with Jacobian at the left hand side) without referring to the above symplectic techniques. 
They successfully use it in Nambu formulation of a class of systems whose symmetry algebras close 
into some simple Lie algebras (see the next section) and then they develop, in analogy to the 
classical case, interesting quantum versions of NB. 



7 



the first two pairs arc fixed), and so on. In aU this cases the signs of permutations are the 
same since each is obtained from the initial one by interchanging two fixed pairs of indices. 
Hence, in the right hand side of (23) there are n\ identical copies of each non-zero term. 
Thus, (23) can be written, in view of identifications given by (24) , as 

^""i^h, $f,J=n\K{n){f,H}, (25) 
where K{n) represents all different partitions which are non-zero 

K{n) = Ne'^ - ■ ■ ■ A„_J. (26) 

Since e^-'^~^ — 1, the factor N is found to be (— l)("-i)("-2)/2 -^y. computing the parity of 
permutation 

/123 4 5 6 7... 2n-l 2n\ 
\^123n + 2 4n + 35... n + 1 2n J ' 

As a result, by defining B^j — {Hi,Aj} and the (n — 1) x (n — 1) matrix B — (By) with 
determinant 

detB^e'' Bu, ... 

Eq. (25)can be written as 

• • • , €/.J = (-l)("-^)("-^)/^n!cietE{/, H}. (27) 

It will be convenient to give the explicit calculation of (27) for n = 1, 2, 3. For n = 1 
we immediately get ^{^f^, ^/j) = {f,H}, and detB — 1. In the case of n = 2 we obtain 
directly from (23) 

J^'(l/l, l/4) = (-irH/l,/2}{/3,/4} + (-l)^H/3,/4}{/l,/2}, 

where ei is the parity of identity permutation and 62 is the parity of 

(3 4?2)=(13)(24) ^ 62 = = 61. 
Hence 0^(^/i, . . . , ^/4) = 2{/, HjBn, and detB = Bu = {-^"1, ^1}. Forn = 3 we have 

^'i^f., Oa)={/l,/2}E{[(-l)^'^ + (-l)^?]{/3,/5}{/4,/6} + 

i=l 

where €{ ;j — 1, ,2, 3, 4, are for the following permutations 



1 2 3 4 5 6 

1 2 3 5 4 6 

1 2 3 4 5 6 

1 2 4 6 3 5 

1 2 3 4 5 6 

1 2 3 6 4 5 

1 2 3 4 5 6 

1 2 4 5 3 6 



(45) ^ 6} = 1 , 
(3465) ^ 6^ = 1 , 
(465) ^ = , 
(345) ^ = . 
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For 62 and e|, the first pair (12) niTist be interchanged, respectively, with the second pair 
and the last pair in the above permutations; in both cases we have the same signs as e\. 
Hence Eq. (28) can be written as 

^'{^h, ■ ■ ■ , ^/s) = -3!{/, H}{B^,B22 - B,2B2^), (29) 

and we get detB = B11B22 — B12B21. 

In view of Eq. (18) a comparison of Eqs. (22) and (27) yields 

d{f, H, Hi, ...Hn-i, Ai, An-i) , 



In terms of 



(-l)"+MeiS{/,//}. (30) 



d^'x' ^dx^ N ■■■ A dx'-^ A A • • • A dx" , 

we obtain for F, by comparing (11) and (30) 

n 

F = detB YjK-^T\Qt^d''q A - dp. Hd''q' A dy], (31) 

1=1 

where the summation is explicitly written to avoid any confusion. Note that *{df A F) is 
equal to the left hand side of Eq. (30). We conclude this section by the following statement. 

li H ^ constant then the constants of motion H, Hi, Ai & Slni'^n — 1) are functionally 
dependent (F = 0) if and only if detB — 0. These functions are functionally independent 
(F ^ 0) where detB ^ 0. 



V. NAMBU FORMULATION OF MAXIMALLY SUPERINTEGRABLE 
SYSTEMS AND THEIR SYMMETRY ALGEBRAS 

Since the Jacobian determinant (30) is proportional to {/, H}, the correct time evolution 
of / can be expressed by the properly normalized NB as in Eq. (15) with the normalization 
coefficient C = {—lY^^ /detB. Indeed, the bracket 

{-lY+^d{f, H, Hi, Hn-l, Ai, An-l) 



detB d{q^, g", pi, p„) 

written in terms of (2n — l)-form 



detB 

produces the correct Hamiltonian time evolution 



(32) 



r' - ^i^r (33) 



f^-{LH}P^ = {f,H}, (34) 
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where / is an arbitrary function. Therefore, every nD maximal superintegrable time- 
independent Hamiltonian system defined by (1) admits equivalent Nambu formulation. 

It must be emphasized that as we have proved in Sec. Ill the normalized NB defined, 
for two arbitrary functions h, f E A, hj 

{f,h}^^^^=*idfAdhATH), (35) 
obeys the Jacobi identity. Here the (2n — 2)-form Th is defined as 



H 



detB 



-dHi A ■ ■ ■ A dHn-i A c/Ai A ■ ■ ■ A dA^-i. (36) 



The above discussion singles out an important special case in which the symmetry alge- 
bra of a 5'/f/(2n — 1) system is such that detB is everywhere a non-zero constant. In such 
a case the constants of motion are globally functional independent and the Nambu formu- 
lation is possible without any (nontrivial) normalization coefficient. For this reason the rest 
of this section is devoted to a discussion of some general points of symmetry algebras of 
superintegrable systems. 

Besides the vanishing PBs given by (1) the defining relations of the symmetry algebra of 
a SIhC^ti — 1) system contain the following 

{Hi, Aj} = Bij, {Ai, Aj} = Cij. (37) 

By the Jacobi identity Bij and Cij are constants of motion and each is functionally dependent 
to the constants of motion specified by the set Slui'^n — 1). Hence, each of them can be 
expressed as X = X(iJ, Hi, Aj). By making use of the identity 

re-l 

{R, X} = {R, H}dHX + ^ {{R, H.i}dH,X + {R, A,}dA,X) (38) 

i=l 

we obtain 

n-1 

{H,,X} = Y.B,,dA,X, (39a) 
i=i 

n-1 

{A,, X} = Y. {-B,,dH,X + C,,dA,x) , (39b) 
i=i 

Obviously, the symmetry algebra of Slni'^n — 1) system is a Lie algebra if and only if each 
of Bij and Cij is at most first order in H, Hi and Aj. If each of Bij and Cij is a polynomial 
of degree at most > 1, then the symmetry algebra is called to be a polynomial Poisson 
algebra of degree k. There may also be cases in which the right hand sides of Eqs. (39) are 
polynomials but some of Bij and Cij are not. In such cases, by including the non-polynomial 
ones into the set of symmetry algebra we again obtain polynomial algebras. Although these 
possible cases are by no means exhaustive, especially for low values of n they are likely to 
occur as interesting structures [|T9| , ^ , ^ , ^ . 
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VI. MULTI-HAMILTONIAN STRUCTURES OF MAXIMALLY 
SUPERINTEGRABLE SYSTEMS 



To ease the calculation of this section we shall use the notion Hk] /c = 0, 1, . . . , 2n — 2 
such that Ho — H and Hn+i-i — Ai, i — 1, . . . , n — 1. Then, in terms of {2n — 2)-forms 

= JetB A A • • • A dHk-i A dHk A dHk+i A • • • A d//2n-2, (40) 

where a hat over a quantity indicates that it should be omitted, we can define 2n — 1 different 
normalized NBs as follows 

{f,h}^Pp^* {dfAdhATH,), A; = 0, 1, 2n-2. (41) 

Each of these brackets gives the original time evolution provided that we choose the new 
Hamiltonian function to be H]^: 

f^-{LH,}^N^P = {LH}. (42) 

In such a case, the system given by (1) has the so-called multi-Hamiltonian structures 
property: it can equally well be described by any one of the 2n — 1 pairs (iJ^, {, j^vp)- 

We shall now prove that the above defined brackets are pairwise compatible, that is, 
{,}n^p^^ — Itvp + b{,}^^p satisfies the Jacobi identity for all a, 6 e R, independently 
from the form of detB. Evidently, there is no loss of generality in taking ki — Q,k2 = 1 and 
a = 1 = 6. Then let us consider 

{/, {9. h}f^}f^, +cp^ [{/, {g, h}Pp}^^^p + cp] + [{/, {g, + cp], 

= {*[df A d*{dg AdhA Tho) A Th,] + cp} + (43) 
{*[df A d*{dg AdhArH,)A Th,] + cp}, 

for three functions /, g and h. In writing Eq. (43) we have made use of the fact that {, }j^p 
and {, satisfy the Jacobi identity separately. The definition of {, }j^p requires 2n — 1 
independent functions but (43) involves 2n + 2 functions. Therefore, in the most general 
case any two of (/, g, h), say / and g, must depend on the other 2n independent functions. 
That is, we can take 

2n-2 2n-2 

df = aidh + J2 hdHk, dg = a2dh + ^ CkdHk, 

k=0 k=0 

where ai, 02, bk, Ck are arbitrary constants. When these are substituted in (43), with special 
care paid to the signs, the right hand side vanishes. 

To point out another important property of superintegrable systems we shall use local 
canonical coordinates x''; 7 = 1, . . . ,2n such that = q\ x""*"^ = pj for j — 1, . . . , n. 
Then, in terms of the so-called Poisson tensor components 

^ap ^ t^r^^ap^o ...f....T.„-.^^ . . . dJh-ldlhdHj^ dH2n-2 

^''^ detB dx'^o ^xt^-i dx^^k dx'^k+i dx^^^-i ' ^ 
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we rewrite Eq. (41) as 



{/'MSJ!. = A^i5^|^, k = 0,l,...,2n-2. (45) 

In Eqs. (44-45) and below summations over repeated Greek letters range from 1 to 2n. It 
is not hard to verify that in terms of Poisson tensor the Jacobi identity means that 



^iii^^i') + ^P = 0^ (46) 



where cp indicates the cyclic sum with respect to superscripts r], a and p. We shall now 
prove that all of these Poisson tensors are singular, that is, A(fc) being the 2n x 2n matrix 
with elements A°^^ we have 

ciet(A(,)) = A;,7 . . . A?-- = 0. (47) 

The easiest way to prove Eq. (47) may be as follows. Let {ea '■ a = 1, . . . , 2n} be a 
basis of vector space R^" and let us define 2n vectors v*^") = A^g'^e^, where , without any 
loss of generality, we have taken k = 0. Then for k' = 1, . . . , 2n — 2 we have 

because of the contraction of two symmetric and two antisymmetric indices. This proves 
that each of v'-"-* is perpendicular to the set of 2n — 2 linearly independent vectors Sh = 
{'VHk' : k' = 1, . . . , 2n — 2}. Hence the 2n vectors v*^"^ are linearly dependent and 
therefore the matrix of their components, which is the matrix A(o), is singular. In fact the 
rank of this matrix is two since this is the dimension of the orthogonal complement of the 
set Sh- 

Evidently, all of the Poisson tensors defined above are pairwise compatible, but as they 
are singular they do not lead to any symplectic structure. 



VII. APPLICATIONS 

A. Calogero-Moser system 

The Calogero-Moser system is one of the four nD systems which are known to be max- 
imally superintegrable for any finite integer n |^,|1^,[T6|. The other three systems are the 



Kepler-Coulomb problem, harmonic oscillator with rational frequency ratios, and Winter- 
nitz system. Here we shall consider n = 2 (two particles) case of the Calogero-Moser system 
described by 

where g is a, constant. Constants of motion for Hcm can be written as 
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m 

A, = (49) 

m 

They obey the following relations 

{HcM.Hi} = {HcM,M} = {HcMi^'i} = 0, 

{A,,A[} = 0, {H,,A,} = 2A[, {H,,A[} = -2A,A[. (50) 

The relation {Ai, A[} — imphes (and is implied by) the fact that Hcm, Ai and A[ are 
functionally dependent. From Eqs. (48) and (49) one can easily identify this dependence as 
A[ = {AHcM-rnAl)/2m. Thus, functional independent set we can take {Hcm, Hi, Ai), 
or {Hcm, Hi, A[). In the former case the symmetry algebra is spanned by {Hcm, Hi, Ai) 
and since 

Bn = {H,, A,} = -Al + -Hcm, (51) 

m 

it is a quadratic Poisson algebra. In the latter case the symmetry algebra is spanned by 
{Hcm, Hi, A[, B{-^), where = {H^, A'^} such that 

B[\ = -8Af + -A['HcM. (52) 
m 

From Eqs. (39a) and (39b) it is found that 

^^^^ = 2 ^ = + ^^^^^ 

K,i?ii} = -^^ = 0. (53b) 

Hence we have again a quadratic Poisson algebra. 

In both cases the time evolutions can be written in the Nambu and Hamiltonian me- 
chanics equivalently as 

"^■^ ^ {f,HcM,H,,A,} (54a) 



dt Bu 
1 



B 



-{/, Hcm, H,,A[} = {/, Hcm} ■ (54b) 



Both of the brackets 
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{f,9}Pp = -^{f,g,H„A,}, 

-Dll 

{f,9}%%^-^{f,9,H,,A[}, 
Hu 
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defined for two arbitrary functions /, g satisfy the Jacobi identity. This is also the case for 
the normahzed NBs corresponding to the Hamiltonian functions Hi 

{/> 5}Svp = -E~{f^ 9, HcM, A} , 

-Dll 

{f,9}%i = ^{f:9:HcM:A[}, 

and to Ai and to A[ 

{f, 9}np = -^{/' 9, HcM, Hi} , 

-Dll 

{/> 9}'np = ~-Br{f^ 9, HcM, H'l} . 

-Dll 

Finally in this subsection we should note that by redefinition of the constants of motion 
we can make the normalization coefficients trivial at the expense of restricting their domains 
of definition. As an example we consider the constant of motion 



2h Ai + h' V m ' 

which satisfies {Hi,A'l} = {Hi, Ai}dAj^A" = —1. Then, we can rewrite Eq. (54a) as 

I = {/, HcM, Hi, A"} = {/, Hcm}. (55) 
Two more alternative brackets can be defined as 

{f,9}f^ = -{f,9,HcM,A';}, 

{f,9}T^={f,9,HcM:Hi}. 

Like others, all these brackets obey the Jacobi identity are pairwise compatible and each 
is degenerate. Explicit forms of corresponding Poisson tensors can also be found as will be 
done for the next application. 



B. Landau problem 

We now consider the well-known Landau Hamiltonian if^, and first establish its super- 
intcgrability and Nambu formulation. Explicit expressions of three different Hamiltonian 
structures of this problem will be presented in the next subsection. 

For a particle of charge q > and mass m moving on the g^g^-plane under the influence 
of the perpendicular static and uniform magnetic field B = dqia2 — dqiai , H^ is (in the 
Gaussian units) 

l^f = \m{vl + vl), (56) 
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where c is the speed of hght, a = (oi, 02) is the vector potential and v = (p — -a)/m is the 
velocity vector [|^. Components of v obey {^1,^2} = qE/rn?c for any B. 

When B is constant the most general form of the vector potential is a = (i3/2)(— g^, q^) + 
VgX ) where x = is an arbitrary gauge function. In such a case we have, in any gauge 
X, two constants of motion 

Hi = m{v2 + ojq^), Ai = -m{vi - ujq^), (57) 

where uj = qB/mc is the cyclotron frequency, (ifi, Ai)/muj correspond to coordinates 
of the cyclotron centre and they satisfy the gauge- independent relations {vj.Hi} = = 
{vj.Ai}, j = 1,2, and 

{Hl, Hi} = = {Hl, Ai}, {Hi, A^} = -mu. (58) 

These relations explicitly show that Hl, Hi and Ai are functional independent constants of 
motion and they close into a Lie algebra structure which can be identified as the centrally 
extended Heisenberg-Weyl algebra. This completes the superintegrability of the Landau 
problem. 

Let us now consider the 3-form F = {muj)~^dHL A dHi A dAi, which can be written as 

r = {vidq^ — V2dq^) A [dpi A dp2 {dpi A da2 — dp2 A dai)] + mu\ ■ dp A dq^ A dq^. (59) 

Making use of this we immediately have 

*((iqAr) = v, 

*{dp AT) = -{viVai + V2Va2) . 

c 

It is straightforward to check that the right hand sides of these equations are the right 
hand sides of the canonical Hamiltonian equations for Hl. Therefore, we can write them 
collectively as 

^u = *(duAr), (60) 
where u = {q^, q^, pi, P2). One may also write Eq. (60) for any function /. 



C. Three Hamiltonian Structures of Landau Problem 



Since Bu = {Hi, Ai} = —muj we can define the following three different Poisson tensors 

^af3 ^ ^ ^a/37172 OHi 

muj 9x^1 dx'y^ ' 

for the Landau problem. These are all computed and in terms of 2 x 2 matrices 



15 



and 

/ = detZ , li = vidiai + f 2^102 , h = Vid2ai + V2d2a2 , (63) 
they can be expressed by the following 4x4 matrices 

1 ^ ] 

B\-Z HY r 



A(0) = -'o + 



A(i) = -V2A(o) + I % ) , (64) 



A(2) = t^iA(o) + 



^2 

-^2 



Here Z denotes the transposition of the matrix Z and in terms of 2 x 2 zero matrix and 
unit matrix 1 

^°=(-lj)' 

stands for the standard 4x4 symplectic matrix. 

We first should note that these Hamiltonian structures are valid in any gauge x- Then, 
as applications of general statements proved in the main text one can explicitly verify the 
following. 

(i) Each of A(fc) provides the correct Hamiltonian equations for the Landau problem which 
can be cast in the following matrix forms 

ju = A(fc)Vi/fc ; = 0, 1, 2 , (66) 

where Hq = H^, H2 = Ai and Viffc is the column matrix of gradient H^. 

(ii) detK{k) = and in fact each of A(^k) has rank two. 

(iii) All A(/(.) satisfy the Jacobi identity (46) and they are pair-wise compatible. 



D. Smorodinsky-Winternitz potentials 

We now consider a set of four potentials first found by Winternitz and co-workers 



Relaying on the assumptions; (i) Hamiltonians are of potential form, (ii) integrals of motion 
are at most quadratic in momenta, they found the following four potentials 



yd) 



= coiAql + ql) + a2qi + ^ , (67b) 

V^'' = ^U + ^ + ^]^ (67c) 
Zr \ r + qi r — qi J 

= _L(c, + «4V^qr^ + /?4yf^), (67d) 
2r 
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where k,Lj,K,a and are some real constants and r'^ = qf + g|. All these potentials 

accept separation of variables in at least two coordinate systems and each has, in addition 
to Hamiltonian function H^^^ = {p^/2m) + V^^\ two constants of motion Hi\Ai^;j = 
1,2,3,4. That is, all the Smorodinsky-Winternitz systems are superintegrable and they 
contain two-dimensional harmonic oscillator and Kepler-Coulomb problem as special cases. 
D-dimensional version of (67a) for D > 2 is known as the Winternitz system HTSj . 



Explicit expressions of constants of motion and their nonvanishing PBs are given alto- 
gether in the Table I. In addition to that given in third column of the table one must add 
the following relations 



(68) 



to the defining relations of symmetry algebras. The constants of motion b[{^ = {Hf' , A^'} 
are found to be 



r(i) aU)' 



B 



(1) 
11 



4 
m 



Bii = -—q2P2 {8ujqi + ^2) - — ( — - 2u;g2 + 



h kqiq2 - ai — - + fJ 

\ m / qi^ 

2pi (pI 



q2' 



W2 



m 



m \m 



5(3) 
-Oil 



2pi^^ ^ 2g2^ fi^^ «3 ^ Ps 



^2 

+ — qp 

mr 



(69) 



"3 



mr \2 r + qi r — qi 
where B[f is not written due to its length and we have defined 



(r + gi, 



- aj ±(3j, L = p2qi - Piq2- 



(70) 



Since s are cubic polynomials in the momenta they can not be written as polynomials 
in the constants of motion [^. But their squares can be expressed as follows 



o(l)2 
-Oil 
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m 



^^(//W 2 _ 2/7«/7(i) + kA^^-) _ 27+ A;) + 4«iif(i 



)2 



16 



^7r(2if«if«-7rA:), 



m 



B 



(2)2 
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r(3)2 
-Oil 



B 



(4)2 
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m 
1 

m 



|(2)^o-(2) 



(2)2 



r(2) 



1 

m 



■^(3) ^4^(3) 2 _ 8^(3)^(3) + 8^+^(3) _ ^2) _ (nA^^ + ^~ H^^^) + K^^t 

Hi'\2H^'^Hi'^ - l7l74-) + A^^\2H('^Af - a,(3,) - 



As a result, the symmetry algebras for all the Smorodinsky-Winternitz potentials are five 
dimensional quadratic Poisson algebras generated by H^^\ Hi\ Ai \ Bjj and 1. Finally 
in this section we will be content with writing the normalized NB of this system as 



-On 



(71) 



It is also straightforward to write out the alternative Poisson structures for each member of 
this class of potentials. 
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TABLES 



TABLE I. Constants of motion and nonvanishing PBs for the symmetry algebras of the 

Sniorodins]<;y-A\^intornitz ]:)()tontials. The alibrcn'iations "J and L -axc dofinod b\- Eq. (70). 





Constants of Motion 


Nonvanishing PBs 




H? = 'i + H + ^^ 


{i7f),5S)} = -A[F«(,fW-2ilW) 
+2^(^(1) -7i+)] 




H?^ = ^p? + Awql; + 0201 


{Af\Bf^} = -^(3ijP' - 4iJP)ijP + ij(2)^ 


ij(3) 




{4^5i?} = ^[4^^'-2i^(^)(2i?P-73^) 


i^{4) 


+^i>2 + §(r + giVr-9i 

= («4V^ - gi - ^A^r + qi) 





20 



